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Abstract. We have used VLT/UVES to spatially resolve the gas disk of /3 Pictoris. 88 extended emission lines are observed, 
with the brightest coming from Fel, Nal and Call. The extent of the gas disk is much larger than previously anticipated; 
we trace Nal radially from 13 AU out to 323 AU and Call to heights of 77 AU above the disk plane, both to the limits of 
our observations. The degree of flaring is significantly larger for the gas disk than the dust disk. A strong NE/SW brightness 
asymmetry is observed, with the SW emission being abruptly truncated at 150-200 AU. The inner gas disk is tilted about 5° 
with respect to the outer disk, similar to the appearance of the disk in light scattered from dust. We show that most, perhaps all, 
of the Nal column density seen in the 'stable' component of absorption, comes from the extended disk. Finally, we discuss the 
effects of radiation pressure in the extended gas disk and show that the assumption of hydrogen, in whatever form, as a braking 
agent is inconsistent with observations. 
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1. Introduction 

The young, near-by main-sequence star /3 Pictoris has been the 
subject of intense studies ev er since it was discovered to har- 
bour circumstellar cold dust ( Aumann fun, distributed along 
a linear shape ( Smith & Terrile 1984), interpreted as a "de- 
1 bris disk" jBackmari&ParesceTl 993). These studies have been 
largely motivated by the possibility of observing an analogue to 
the solar system in its early stages, in the hope of finding clues 
to the mechanisms of planet formation. Asymmet ries found 
in the disk from light s cattered by the dust jKalas & Jewittl 
1 19951 iHeap et all l2.000h have indeed been suggested to be 
the signature of per turbing planet(s) (Mou llietetaT]ll997l 
lAueereau et alJ 1200 lh . but so far no direct detection of a 
planet around /3Pic has been made. Asymmetri es have also 
been d etected in thermal emission from the dust jLiseau et alJ 
120031 IWeinberger. Becklin. & Zuckermanl 120031 with refer- 
enc es therein). R e cent pa p ers reviewing the P Pic dis k are those 
by lArtvmowicj J200C ). lLa grange. Backman & Artvmowicj 
(2000) and Zuc kermadfeOOll) . 

Circumstellar gas, see n in absorption against the star, 
was also found early on ( Hob bs et alJ fl98 5). thanks to the 
favourable edge-on orientation of the disk. Finding and char- 
acterising the gas content is important for understanding its 
relation to the dust and the general evolution of the disk 
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dArtvmowic jl200oT) . Gas is also useful as a probe of physical 
conditions in the disk, where density, composition, temperature 
and bulk velocities under favourable conditions can be directly 
estimated. 

The gas found at relative rest to jSPic, consisting of met- 
als, raised the problem why it is not blown away from the 
system by the high radiation pressure. Lagran ge et alJ {1998) 
made some detailed calculations and found that the gas drag 
from a dense enough HI ring (rim ^ 10 5 crrT 3 ) close to the 
star (~0.5 AU) could brake migrating particles sufficiently, pro- 
vided they started out inside the ring. The picture was com- 
plicated by the announcement of H? detected in emission 
by the Infrared Space Observatory dThi et al.ll200lh . implying 
large quantities (~50M e ) of molecular hydrogen, and the sub- 
sequent report of sensitive upper limits (A^Pb) 55 10 18 crrT 2 ) 
of H2 from F UV absorption lines using /?Pic as a back- 
ground source (Lecavelier des Etangs et al. Il200ll) . In addition, 
lOlofsson. Liseau. ^^randeke^l200 l[ hereafter Paper I) found 
spatially resolved widespread gas emission from Nal in the 
disk, stretching out to at least 140 AU. 

Here we present observations improved by a factor of two 
in both spatial and spectral resolution, as well as a greatly in- 
creased spectral coverage, compared to Paper I. We put empha- 
sis on the observed spatial structure of the gas disk, derive an 
empirical density profile of Na I atoms and use a photoionisa- 
tion code to construct disk models consistent with our observa- 
tions. We discuss implications of the radiation pressure under 
various conditions derived from these models. Results from a 
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detailed study of the observed chemical abundances will be dis- 
cussed in a forthcoming paper. 

2. Observations 

The star ft Pictoris was observed with the echelle spectro- 
graph UVES on the 8.2 m Kueyen telescope at the Very Large 
Telescope of ESO, Paranal, Chile. The observations were made 
in service mode to take advantage of the excellent seeing con- 
ditions occasionally provided by the site. Service mode means 
that the observers prepare observation blocks that are executed 
by the local telescope operators provided certain conditions 
(like air mass and seeing) are met. A total of 12 observation 
blocks during 2001-2002 (see Tabled lasting about an hour 
each were successfully executed. The air mass and seeing at 
zenith were always below 1.5 and 076, respectively. The UVES 
spectrograph was used in standard dichroic modes, meaning 
that a dichroic mirror was used to split the light into two wave- 
length ranges, each using its own echelle grating and CCD de- 
tectors. Two types of gratings were used for each arm, together 
covering the wavelengths 3 300 A to 10400 A with a spectral 
resolution ranging from /?~80000 in the blue to R~l 10000 
in the red. The blue arm was equipped with a single CCD of 
type EEV 44-82 and size 4096x4096 pixels, while the red 
arm had a mosaic of two CCDs of sizes 2 048x4096 each, of 
the types EEV 44-82 and MIT-LL CCID-20. The spectroscopic 
slits were 074 and 073 wide in the blue and red arm, respec- 
tively. To adequately separate the echelle orders, slits were re- 
stricted to lengths between 8" and 12", depending on setting 
(see Table 0. Slits were positioned at 4 overlapping locations 
along the ft Pic dust disk, and at 4 locations orthogonal to the 
disk, reaching a distance of 17" from the star in the disk plane, 
and a height of 6" above the plane (see Fig.Q. From the data, 
we determined the positioning of the slits, with respect to ft Pic, 
to be better than ~072. Calibration data (fiat-fields, bias frames, 
Th-Ar lamp wavelength calibration spectra, flux standard stars, 
etc.) were provided by the UVES calibration plan. 

3. Data reduction 

3.1. Pipeline reduction 

The data were reduced using a modified version of the UVES 
pipeline 1.2, that runs in the ESO MIDAS environment. The 
pipeline automatically generates the appropriate calibration 
sets from the calibration data obtained as a part of the UVES 
calibration plan, and apply them to the science frames. The raw 
frames are bias subtracted, flat-fielded, background subtracted, 
order extracted, and wavelength calibrated. An error was cor- 
rected in the pipeline version 1 .2, that associated a wavelength 
calibration to the wrong order when 2D extraction was selected 
without merging the orders. 

The orders of the echelle grating are slightly inclined on the 
CCD with respect to pixel columns. Since the observations with 
the star on slit were performed under excellent seeing (the see- 
ing measured in the spectra was ~077), this resulted in a very 
high flux gradient between the pixels on and off the star. By de- 
fault, the orders are resampled linearly, but this caused the spec- 
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Fig. 1. Orientations and positions of the 8 different slit settings 
used. The background image shows obser vations of the dust as 
obtained by HST/STIS (Hea oet all2o"o(il) . 



trum to show a periodic 'spiky' pattern in the flux. We found 
that linear interpolation does not work properly due to the 
strong non-linear shape of the point spread function (PSF) at 
the pixel resolution. By modifying the order extracting routine 
to use a third order Catmull-Rom spline interpolation instead 
of linear, we found the spikes to be greatly reduced (though not 
completely eliminated). 

The reduction procedure was complicated by the sheer load 
of data. With about 600 spectra, each 4096x4096 pixels, plus 
calibration frames etc., the raw data were close to 50 GB. 

The end product of the pipeline reduction consisted of in- 
dividual orders separated into individual files with the spatial 
information along the slit preserved, and a wavelength calibra- 
tion guaranteed to be better than 0.5 km s _1 . 

3.2. PSF subtraction 

To trace the circumstellar gas emission as close as possible to 
the star, we had to subtract the scattered light from the star, in 
the cases where the star was in the slit. To do this, we estimated 
the stellar spectrum S by centering a small aperture on the star. 
To allow for slow gradients, we took advantage of the fact that 
the gas emission lines are narrow (a few pixels) and produced 
a median filtered version M of the stellar spectrum S , using a 
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Table 1. VLT/UVES observation log 



UT Date 


/?Pic a 


Grating b 


Width c 


Length 


Offset d 


Orientation d 


Exposures 


2001-12-02 T 03:38 


ON 


437 


074 


10" 


-474 


Parallel 


69 x 15 s 


2001-12-02T03:38 


ON 


860 


073 


12" 


-474 


Parallel 


47 x 25 s 


2001-12-02T05:31 


ON 


390 


074 


8" 


-374 


Parallel 


62 x 20 s 


2001-12-02T05:31 


ON 


580 


073 


12" 


-374 


Parallel 


48 x 25 s 


2002-01-03 T 06:21 


ON 


390 


074 


8" 


-374 


Parallel 


62 x 20 s 


2002-01-03 T 06:21 


ON 


580 


073 


12" 


-374 


Parallel 


48 x 25 s 


2002-01-06 T 05:51 


OFF 


390 


074 


8" 


-370 


Orthogonal 


3 x 1205 s 


2002-01-06 T 05:51 


OFF 


580 


a/ 3 


12" 


-370 


Orthogonal 


3 x 1195 s 


2002-01-16 T 02:47 


OFF 


390 


074 


8" 


-11" 


Parallel 


3 x 1205 s 


2002-01-16 T 02:47 


OFF 


580 


a' 3 


12" 


-11" 


Parallel 


3 x 1195 s 


2002-01-16 T 03:51 


OFF 


390 


074 


8" 


-670 


Orthogonal 


3 x 1205 s 


2002-01-16 T 03:51 


OFF 


580 


073 


12" 


-670 


Orthogonal 


3 x 1195 s 


2002-01-19 T 03: 11 


ON 


390 


074 


8" 


374 


Parallel 


62 x 20 s 


2002-01-19 T 03: 11 


ON 


580 


a/ 3 


12" 


374 


Parallel 


48 x 25 s 


2002-01-19 T 04:25 


ON 


390 


074 


8" 


374 


Parallel 


62 x 20 s 


2002-01-19 T 04:25 


ON 


580 


a' 3 


12" 


374 


Parallel 


48 x 25 s 


2002-02-09 T 03:55 


OFF 


390 


074 


8" 


370 


Orthogonal 


3 x 1205 s 


2002-02-09 T 03:55 


OFF 


580 


a/ 3 


12" 


370 


Orthogonal 


3 x 1195 s 


2002-02-10 T 03: 19 


OFF 


390 


a/4 


8" 


11" 


Parallel 


3 x 1205 s 


2002-02-10 T 03: 19 


OFF 


580 


073 


12" 


11" 


Parallel 


3 x 1195 s 


2002-02-11 T 04:04 


OFF 


390 


074 


8" 


670 


Orthogonal 


3 x 1205 s 


2002-02-11 T 04:04 


OFF 


580 


a/ 3 


12" 


670 


Orthogonal 


3 x 1195 s 


2002-09-24 T 08:09 


OFF 


437 


a/4 


10" 


-370 


Orthogonal 


3 x 1205 s 


2002-09-24 T 08:09 


OFF 


860 


a/ 3 


12" 


-370 


Orthogonal 


3 x 1195 s 



a ON means that the star was inside the spectroscopic slit, OFF that it was outside. 

b The gratings correspond to the following wavelength ranges: 326^145 nm for grating 390, 373-499 nm for grating 437, 476-684 nm for 

grating 580, and 660-1060 nm for grating 860. 
c Width and Length refer to the dimensions of the slit. 

d Offset is the centre of the slit relative to /3 Pic in the plane of the disk, where positive offsets are to the north-east. The orientation of the slit 
is relative to the disk plane. 



window of 30-60 pixels. We also obtained the median filtered 
spectrum m, at each spatial position ;, using the same median 
filter window. The stellar PSF subtracted spectrum was then 
estimated as Si = r,- - m,-5 /M, where r; is the spectrum at po- 
sition i. For the observations where the star was off the slit, we 
just subtracted the median filtered background scattered light 
(scattered both locally by the dust disk and in the atmosphere / 
telescope), s, = r, - m,-. 

3.3. Flux calibration 

To flux-calibrate the data, we made use of the master re- 
sponse curves provided by the UVES team. The master re- 
sponse curves are generated from long-time monitoring of the 
sensitivity trends of the instrument, and are provided for vari- 
ous standard settings and periods in time. The claimed absolute 
flux calibration is 10 %, but comparing with standard flux stars 
observed (with an open slit) the same nights as (3 Pic we found 
the derived fluxes to deviate by as much as 40 % from tabu- 
lated values for these stars. We thus used the master response 
curves to correct only for the over-all sensitivity dependence 
on wavelength, and used the standard star spectra obtained in 
connection with the /3 Pic observations to correct the absolute 
flux calibration. By looking at the variations of the derived ab- 



solute fluxes of disk emission lines observed several times, we 
estimate that the error on the absolute flux is about 5-10%. 
Since we use a narrow slit (0'.'3-0'.'4), the slit losses from a 
point source like [} Pic are substantial and very seeing depen- 
dent. The error of the absolute flux from the star is therefore 
expected to be substantially higher. Our main concern, how- 
ever, is the flux of the disk emission, which, due to its spatial 
extension, is much more stable with regard to seeing variations. 

3.4. Heliocentric wavelength correction 

The wavelength calibration obtained from the UVES pipeline 
does not correct for the velocity of the instrument relative to 
the centre of the Sun. To transform the observed spectrum to a 
heliocentric frame, we made use of the software program RV 
written by P. T. Wallace and C. A. Clayton 1 . The obtained ac- 
curacy in the transformation is better than 0.01 km s~', fully 
adequate for our purposes. 

3.5. Line flux measurements 

From the reduced, PSF subtracted and wavelength calibrated 
spectra, the disk emission lines were measured in several steps. 

1 RV is available from http://star-www.rl.ac.uk 
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First the lines were identified by correlating their measured 
wavelengths with the atomic line database provided by NIST 2 . 
This procedure was greatly simplified by the high accuracy of 
the wavelength calibration, better than 0.5 km s -1 , correspond- 
ing to 0.01 A at 6 000 A. A median systematic velocity was then 
calculated for the brightest and most accurately measured emis- 
sion lines, assuming all shared the same systematic radial ve- 
locity. For the brightest lines, the measured flux as a function 
of apertures centred on the systematic velocity was evaluated. 
A large aperture samples more signal, but also more noise. We 
were therefore interested in finding the best balance, in order to 
achieve the highest signal to noise (S IN). The apertures used to 
measure emission lines were consequently chosen as a function 
of line strength and background noise, with smaller apertures 
for fainter lines. We assumed that all lines from a particular ion 
share the same spatio-spectral profile, and scaled the flux mea- 
sured in small apertures of fainter lines with the ratio between 
equally sized and maximum sized apertures placed on the sum 
of several bright lines. In this way we estimated the signal from 
a faint line without integrating up too much noise. The method 
is analogous to methods used in aperture photometry to mea- 
sure star fluxes in, e.g., a CCD image. 

For the slits placed orthogonally to the disk, we summed 
up all flux in the spatial direction, i.e. along the height of the 
disk (Table 13. The quantity thus derived has the unit of flux 
per arcsecond. 



4. Results 

We detected 88 spatially extended emission lines from the /3 Pic 
gas disk, identified as emission from Fe I, Na I, Ca II, Ni I, Ni II, 
Ti I, Ti II, Cr I and Cr II. Table[2]shows a selection of the bright- 
est lines. Following the brightest emission lines (S /N~50) from 
Nal and Fel radially (Figs.|2l&|3|l, we observe a strong asym- 
metry between the north-east (NE) and the south-west (SW) 
parts of the disk, simil ar to the b rightness asymmetry in the 
dust emission dKalas & Jewiti 19951) . but asymmetric to a much 
higher degree. The NE gas emission extends smoothly to the 
limits of our observations (17", corresponding to 330 AU at 
the distance of /3Pic), whereas the SW emission is abruptly 
truncated at 1 50-200 AU. In the inner regions, the SW emis- 
sion dominates over the NE part, in agreement with Fig. 2 of 
Paper! 

Determining the centre of emission from the orthogonal 
profiles, it becomes apparent that the inner part of the disk is 
slightly tilted with respect to the outer parts, in particular on the 
NE side (Fig.|4|i. We estimate this tilt to be 5°±2° (l<x), similar 
to the 4°-5° tilt o bserved by HST/STIS in the inner dust disk 
llHean et all200(l 

The scale height (FWHM) of the gas disk, estimated on 
observations deconvolved with a Gaussian of 0/7 to simulate 
the seeing, i s -20 AU at 3" (5 8 AU), similar to the dust disk 
scale height (Hea p et al.ll2000l) . At 6" (116AU), however, the 
gas disk is significantly thicker, ~30AU compared to ~15 AU 
for the dust (see Fig.|5Jl. 

2 http://physics.nist.gov 




Radial velocity [km/sj 

Fig. 2. Na D 2 (A& = 5889.95 1 A) emission from the Pic disk, 
as seen through four slits parallel to the disk (see Fig.0. The 
vertical axis is along the spatial direction, with positive offsets 
north-east of the star, while the spectral dispersion is along the 
horizontal axis, centred on the /3Pic rest frame of the NaD 2 
line. The inner 1 1 AU, showing mostly residual noise from the 
PSF subtraction, have been masked out. The velocity shift of 
the NE and SW side is due to Keplerian rotation of the disk, 
with the SW rotating towards us (Paper I). The grey scale has 
been scaled non-linearly with the intensity, in order to bring 
out the bright disk structure close to the star as well as the faint 
features at greater distances. The emission in the NE can be 
traced out to the limits of our observations at 323 AU, while 
the emission in the SW ends abruptly at 1 50-200 AU (see also 
Fig. QJ. Close to -25kms _1 , telluric NaD 2 emission is seen 
covering the slit. The relative velocity of the sky emission to 
/3 Pic varies slightly between the mosaiced observations due to 
the orbital motion of the Earth between the epochs. 
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Fig.3. Determined surface intensity as a function of radial distance in NE and SW for NaDi 2 and Fel /13859.91 1,/13820.425. 
The fit of Sect. 15. H is shown for NaID2 as a an unbroken line. Note the sharp decrease in flux in the SW at 150 AU-200 AU. 



Table 2. Selected gas emission lines, 3" SW in disk 



Line 


air 


Flux" 


flux 


low 


hi g h 


Fel 


3820.425 


61.3 


2.5 


0.859 


4.103 


Fel 


3859.911 


111.3 


3.4 


0.000 


3.211 


Nal D 2 


5889.951 


79.2 


4.0 


0.000 


2.104 


NalDi 


5895.924 


42.4 


2.5 


0.000 


2.102 


CallK 


3933.663 


12.0 


1.0 


0.000 


3.151 


CallH 


3968.468 


16.8 


1.2 


0.000 


3.123 



a Wavelength of transition in air, in units of A. 

b Flux in units of 10~' 6 erg s^'cm^arcsecT 1 ; see Sect. 13.51 

c Lower and upper energy levels of transition, in units of eV. 



A very complex height profile is shown by the Ca II H & K 
lines (Figs.|5]&|6}. Especially at 6" distance from the star, the 
emission from the disk midplane is much fainter than the emis- 
sion away from the midplane. In particular, the detected Call 
emission, at 6" SW of the star, keeps increasing to the limits of 
the spectroscopic slit 4" above the disk midplane, meaning that 
there is a significant number of Ca II ions at 77 AU height above 
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Fig. 4. Measured positions of the Na I and Fe I gas emission rel- 
ative to the plane at position angle 30°75. The bars correspond 
to derived scale heights, obtained by deconvolving the ob- 
served height profiles with a Gaussian of 0/7 FWHM simulat- 
ing the seeing. The inclined dash-dotted lines show the empiri- 
cal opening angle of the gas disk, corresponding to H/r — 0.28 
(see Sect. 15. H . The error in positioning of the slits is on the 
order of a' 1. 



the midplane at 116 AU distance from the star. From Fig.[6]it 
is evident that the radial velocities of these Ca II ions are small 
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Table 3. Radial velocities relative ySPic of ions measured in 
absorption 



Ion 


Av a 








l 


Fel 


-0.0 


0.3 


Nal 


-1.2 


0.3 


CalP 


-0.3 


2.1 


Tin 


0.2 


0.8 


Nil 


0.4 


0.4 


Nill 


2.8 


3 


CrII 


2.2 


3 



The positive direction is radially away from us, towards f} Pic. 
Errors are relative only. An additional systematic error of 
0.5kms~' comes from the uncertainty in the system velocity of 
ySPic. 

The radial velocity of Call is measured on the infrared lines 
,17291.47 and ^7323. 89, and not the optically thick Call H & K 
lines. 



- on the order of a few kms -1 . Implications of this fact are 
discussed in Sect. 15.21 

In order to determine the heliocentric velocity of /3 Pic accu- 
rately, we averaged the heliocentric gas velocities as observed 
on each side at the distances 3" and 6" from the star. In this 
way, we can assess the system velocity independently from the 
lines observed in absorption, which are sensitive to possible ra- 
dial velocities of the absorbers (caused by, e.g., radiation pres- 
sure). We determined the heliocentric radial velocity of /JPic 
to be 20.0±0.5 kms -1 , with the error mainly due to the uncer- 
tainty in wavelength calibration. 

In the spectrum of the star, we measured the absorption 
profiles corresponding to observed ground state emission lines. 
Most lines were found to be at, or close to, the systematic ve- 
locity of [i Pic (TableQJ. 

5. Discussion 

We proceed by inverting the observations of the projected disk 
to a radial density law for Na I, and then extend the result by 
calculating the ionisation structure to find an estimate on the 
total Na density. By making assumptions of the abundances we 
construct two models with distinct H densities, test how well 
the models compare with observations of H I and H2, and ex- 
amine how the kinematics of ions in these models are affected 
by radiation pressure. 



5.1. Sodium gas density profile 

From the spatial emission profiles we estimate the density dis- 
tribution of the observed gas by making the following assump- 
tions: 

1. The disk is axisymmetric and the density of the emitting 
medium is well described by 



where r and h are the cylindrical coordinates describing the 
midplane distance and height over the midplane, respec- 
tively, no is a normalising density, ro is the broken power- 
law break distance, a and b are power-law exponents of 
the inner and outer regions respectively, and a is related 
to the scale height H of the disk at midplane distance r as 
H/r = 2 Vhila. 

The inclination of the disk is close to zero, i.e. we see it 
edge-on. 

The emission is optically thin and isotropic. 

The emitting gas traces the full gas population of that ion. 



n(r, h) — no 



in 



2b 



exp 



(1) 



The observed disk is, of course, not axisymmetric but quite 
asymmetric in appearance (Figs. |2J El&EJ- The difference be- 
tween the NE and SW parts of the inner disk is not that dra- 
matic, though, less than a factor of two in the region where 
most gas is located (30-120 AU). Thus the assumption of axi- 
symmetry is probably not too bad. 

We chose the NalDi ^ lines for the inversion because (A) 
the S/N is high (-50), (B) the line ratio between D 2 and D] 
is close to 2, implying that the emission is optically thin, and 
(C) we expect nearly all Na I atoms to be in their ground state 
(see Sect. 15. 3> , meaning the resonance Dj 2 lines trace all of the 
Na I population. The parameters of Eq.^ were fitted by con- 
structing a numerical model, simulating observed spatial pro- 
files as a function of the parameters and ^-minimising the 
difference between the model and the observations. In more 
detail, the model was initiated on a three dimensional grid 
of height x width x depth dimensions 40 AU x 800 AU x 
2000 AU at 1 AU resolution. The central star was assumed to 
emit 4.5 xlO 30 erg s _I A -1 at wavelengths close to the NaID2 
line, and A,-; = 6.22 x 10 7 s~' was used as the D2 transition 
Einstein coefficient for spontaneous emission, in calculating 
the scattering cross-section. For each cell, a specific line lu- 
minosity was calculated, summed up along the line of sight 
and geometrically diluted to the distance of the observer at 
19.3 pc. This 'ideal' image of the disk was then convolved with 
a Gaussian of FWHM = 077, simulating the atmospheric see- 
ing, and sampled along the various spectroscopic slit settings 
of Fig.^to fit the observations. In the fit we used the midplane 
radial emission distribution from 13 AU out to 323 AU of the 
NE, and the orthogonal profiles at 3" and 6" NE of the star. 
The derived parameters of Eq.^for Na I are: 

n = (1.02 ± 0.04) x lO^ctrT 3 
r = 117±3AU 
a = 0.47 ± 0.06 
b = 3.16 ±0.04 
H/r = 0.28 ±0.05, 

where the quoted errors are 1 cr formal fitting errors obtained by 
making an additional number of fits to artificial data generated 
by adding random noise, at the estimated noise level, to the 
observed profile. We have over plotted the Nal fit into Fig. [21 

We can use Eq.^with the above parameters to estimate 
the midplane column density of Nal, to compare with col- 
umn densities estimated from Nal absorption in the stellar 
spectrum seen through the disk. If we integrate only over the 
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Fig. 5. The flux calibrated vertical profile of the disk gas emission at 3" and 6" NE and SW of the star for a few selected lines 
of Na I, Fe I and Ca II. Zero spatial offset refers to the dust disk midplane defined by the position angle 30°75 relative to the star, 
with the positive spatial direction to the north-west (NW). The centre and FWHM for Gaussian fits to the profiles are printed in 
the panels of Gaussian shaped profiles. The SW side dominates the emission at 3", while the NE side slightly dominates at 6". 
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Fig. 6. The velocity distribution of the gas emission as observed by slits orthogonal to the disk at 3" and 6" NE and SW of the 
star for Call, Fel and Nal. To improve the signal to noise, the H & K lines of Call, 10 of the brightest Fel lines, and Dj and 
D2 of Nal were added for each species. The profiles were then normalised so that each level increases one tenth of the peak. 
The noise in the Fe I and Na I plots is below the lowest level curve, while the lowest level curve for Ca II still shows some noise. 
Velocities are referred to the system velocity 20.0km s _1 of /3 Pic, with the positive direction being away from us. Positive offsets 
in the spatial direction are directed to the north-west (NW). The dash-dotted vertical lines correspond to the systematic velocity 
of ySPic, the dotted lines to the Kepler velocity at the projected distance from/? Pic (assuming a stellar mass of M — 1.75 M ), 
while the dashed ellipse in each graph shows the spectral and spatial resolution attained. Note the complexity of the Ca II profile, 
and how the emission at 6" SW (right-most panels) can be traced to the limit of the observations. The detected emission profiles 
from Ni I, Ni II, Ti I, Ti II, Cr I and Cr II cannot be distinguished from Na I or Fe I, due to the lower S/N. 



radii where we are sensitive to Nal emission, that is, from 
13 AU out to 323 AU, we obtain a column density of iV(Nal) = 
3.0 x 10 10 crrT 2 . Extrapolating Eq.^to all radii, that is, from 
zero radius out to infinity, we increase the Na I column density 
by 40 % to JV(Na I) = 4. 1 x 10 10 cirT 2 . This can be compared to 
our observed column density N(Na I) = (3.4 + 0.4)x 10 10 cirT 2 . 
We conclude that most, perhaps all, of the Na I gas seen in ab- 
sorption is situated in the extended gas disk seen in emission. 
We have no reason to doubt that originators of other 'stable' 
gas absorption components also belong to this extended disk. 



5.2. Radiation pressure 

Ions in the /?Pic gas disk are subject to appreciable radiation 
pressure. To evaluate the significance, it is common to de- 
fine the ratio between the forces of gravitation and radiation, 
P = ^rad/^grav Since both the radiation and gravitational fields 
are inversely proportional to the square of the distance to the 
source, /3 is constant throughout the f3 Pic disk. For the observed 
atoms Nal and Fel, ft » 1, meaning that if left alone, grav- 
ity would be largely irrelevant and these atoms would rapidly 
accelerate out of the system at high velocities. This is not ob- 
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served (Tableland Fig.|6j. To investigate what the effects of a 
braking medium are, we solve the equation of motion for a gas 
particle, 



dv _ 
at 



(2) 



where Ff« c is the frictional force. Assuming Fmc = Cv, where 
C is the friction coefficient (proportional to the density of the 
braking medium), and that the distance an atom travels before 
reaching the terminal velocity is much shorter than the size of 
the disk, Eq.[2]may be solved to yield the terminal velocity (as 
t -> oo) 



Vco — ^ ^grav- W) 

Eq.|2 ignores magnetic forces, which should be of little im- 
portance for neutral atoms. If appreciable magnetic fields are 
present, however, ions with net charges may be significantly 
affected. We ignore this complication for the moment and as- 
sume that the magnetic fields are small enough to play a neg- 
ligible role for the kinematics of the observed gas particles in 
the disk. 

For jSPic we estimate /^ai = 250 and jSpei = 18, while C 
depends on the detailed density and temperature structure. We 
calculated terminal velocities as a function of radius for the two 
models of Sects. |5~3~T1 & 15.3.21 Details of the solution to Eq.|5] 
with estimates of C and ft for several ions, are found in Liseau 
<200i. 



5.3. lonisation structure 

To estimate the gas density of all Na atoms in the disk, and also 
get an idea of the total gas density, we need to address the ioni- 
sation structure of the disk. We employed t he one-dime nsional 
photoionisation/PDR code Cloudy ( Ferlan d et al.ll998l) for the 
task. Cloudy consistently maintains the ionization and thermal 
balance, and solves the radiative transfer by making use of the 
Sobolev approxim ation. For a d etailed description of the code 
we used, see lLiseau et al.Ul999h . 

Several assumptions have to be made in order to com- 
pute the ionisation structure, perhaps the most important being 
that of the chemical abundances of the atomic gas. Since we 
have no detailed information, we confine ourselves to study 
two cases of chemical abundances: solar composition and a 
strongly metal depleted model. At first we also considered 
a case with typical interstellar medium abundances, but dis- 
missed it due to the very high degree of Ca depletion compared 
to solar values, in contradiction with observations. We have no 
a priori reason to expect any of these assumed abundances to 
reflect reality, but we believe they still serve as valuable refer- 
ence cases. 

Another important assumption is that of the stellar spectral 
energy distribution (SED). We assume the atmosphere of P Pic 
to be well represented by an ATLAS 9 model dKuruczlll992i) 
with r eff = 8000 K, logg = 4.5 (in cms 4 ) and log(Z/Z Q ) = 
0.0, and the effective radius to be R = l.75R G , implying the 
luminosity L = 11L . Recent VLTI observatio ns of /?Pic esti- 
mate the stellar radius to be (1.735 + 0.128>R o JDi Folco et alJ 



l2003t) . and comparing the ATLAS 9 SED of /3Pic with ob- 
served data from the HST and FUSE archives , we fi nd the 

I — II I 

agreement to be generally excellent (see Liseau 2003). Some 
deviations are found in the far-ultraviolet, where jSPic seems 
to be a slightl y atypical A5 V star with possible chromospheric 
activity jBouret et all2002l) . 



5.3.1. Solar composition disk 

Assuming solar abundances throughout the disk, we can fit the 
radial Na I density of Sect. l5.1l to within 1 % by setting the ra- 
dial density profile of hydrogen nucleii to 



n(H) = 2.25 X 10 3 



x 2.4 / X 5.3 

ro Uo 



(4) 



where ro is the same as in Sect. 15.11 The average fraction 
of neutral Na in the model is 10 3 , while hydrogen is en- 
tirely neutral. The hydrogen column densities are N(HT) = 
8 x 10 18 cirT 2 and N(U 2 ) = 3x 10 I8 ctrT 2 , consistent with 
the observational uppe r limits of A^(HI) £ a few xl0 19 cm~ 2 
dFreudiingetal.il 19951) and N(H ?) $ 3x 10 I8 ctrT 2 (3o% 
iLecav elier des Etangs e t alJl200ll) . The gas temperature of the 
model disk is low, r gas (100AU) = 13 K, while the dust is sig- 
nificantly warmer, rd us t(100 AU) = 50 K. With a radial density 
as in Eq.0]and a vertical profile similar to the Na I gas, the total 
gas disk mass out to 1 000 AU becomes ~0.1 M e . 

Using the equations of motion of Sect. !5.2l to calculate the 
terminal velocity of elements subject to the gas drag of the 
hydrogen density law of Eq.|4]and temperature structure from 
Cloudy, we find that both Fe and Na, both with /? » 1, would 
reach radial velocities on the order of hundreds to thousands 
kms -1 . That is in clear contradiction with the observed radial 
velocities reported in Table [3] The purported explanation of 
Paper I, that Na is ionised most of the time and therefore is 
subject to a small average radiation pressure, does not apply to 
Fe; indeed, both Fel and Fell experience a radiation pressure 
with /3 » 1 . We know of only one remaining plausible expla- 
nation: there is a braking agent keeping the gas from reaching 
high velocities. 

5.3.2. Metal depleted disk 

To investigate the possibility of a more massive disk than impli- 
cated by solar abundances of Nal, we impose a radial density 
profile for hydrogen nucleii, 



n(H) = 10 6 



\2.4 i \5.4 

- At 



(5) 



and tune the metal depletion to reproduce the Nal density 
profile. We found that an overall depletion factor of ~10~ 3 
(more precisely 8.5 x 10~ 4 ) of the metals with respect to so- 
lar composition was required. Due to higher dust-gas interac- 
tion, the temperatures of the gas and the dust are similar over 
the disk, with r gas (100AU) = 36 K and r dust (100AU) = 51 K. 
The column densities of hydrogen are A^(H I) = 6 x 10 20 cm -2 
and N(H2) = 3 x 10 21 cm 2 , significantly higher than the 
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Fig. 7. The terminal velocities for Na I and Fe I in a metal de- 
pleted disk with the radial density of Eq.[5l as a function of 
radius. The horizontal dashed lines show the measured radial 
velocities, as determined from absorption lines (Tabled}, with 
the estimated error shown by the vertical bar. The error is dom- 
inated by the wavelength calibration uncertainty of 0.5 kms . 

observed upper limits reported by iFreudling et al.l (Tl995) and 
iLecavelier des Etangs et al.1 ll200ll) . With a radial density as in 
Eq.[5] and a vertical profile similar to the Na I gas, the total 
gas disk mass out to 1 000 AU becomes ~40 M 9 , similar to the 
mass recently inferred from possible emission in the pure ro- 
tational transitions of H2, J = 2— >0 (28 ami ) and J - 3— > 1 
(17 /im), detected from 13 Pic in ISO/S WS data faii et alJ200ll) . 
Due to the low gas temperature of our model, however, the 
predicted fluxes in the 17/vm and 28 ^m lines are only Fn = 
6 x 10~ 15 erg s" 1 crrT 2 and F28 = 2 x 10~ 15 ergs -1 cm -2 , re- 
spectively, for ortho/para in thermal equilibrium. This is more 
than an order of magnitude lower than the fluxes Fn = 7.7 x 
10 ~ 14 erg s" 1 cm' 2 a nd F28 = 7.0 x 10~ 14 erg s" 1 cm" 2 reported 
bv lThi et"afl d200lh . and can be traced to their higher assumed 
gas temperature of 100K. 

Calculating the terminal velocities of the radiation pressure 
sensitive elements Fe and Na, we produce Fig. [7] where the 
terminal velocity is plotted against the radial distance. We find 
that the derived terminal velocities are consistent with observed 
radial velocities reported in Table|3] 

5.4. Summary 

We have found that in our thin disk, constructed by fitting a so- 
lar composition gas to the Na I observations, the predicted H I 
and H2 column densities are consistent with the observed lim- 
its. With these low densities, however, the radiation pressure is 
predicted to accelerate, e.g., Nal to radial velocities inconsis- 
tent with observations. 



On the other hand, in our thick disk, constructed by fixing 
the H density to be consistent as the braking agent for the met- 
als, the predicted column densities of H I and H2 are both above 
the observed upper limits, although the predicted H2 emission 
is still far below the claimed detection. 

The situation is puzzling, but a possible solution may be 
found if a main braking agent different from hydrogen is in 
action. If, e.g., the observed gas is replenished from dust colli- 
sions rather than being primordial, as seems suggested by the 
close spatial correlation between gas and dust, then a more nat- 
ural candidate might be oxygen or oxygen bearing species (P. 
Artymowicz, private communication). More detailed models 
are in preparation, but hopefully the upcoming Space InfraRed 
Telescope Facility will solve the issue unambiguously by better 
observational constraints. 

6. Conclusions 

Our main observational results are: 

1 . We have observed 88 spatially resolved emission lines com- 
ing from Fe I, Na I, Ca II, Ni I, Ni II, Ti I, Ti II, Cr I and Cr II 
in the f3 Pictoris gas disk. 

2. We trace the gas emission to the limits of our observations, 
from 077 (13 AU) out to 17" (323 AU) radially to the NE 
in Na I, and 4" (77 AU) above the disk plane at radius 6" 
(116 AU) in Call. 

3. The scale height of the gas at 6" from the star is twice as 
high as the equivalent dust scale height. 

4. There is a brightness NE/SW asymmetry in the gas emis- 
sion reminiscent of the dust asymmetry, although much 
stronger. 

5. The inner gas disk is tilted by ~5°, similarly to the dust 
disk. 

6. The heliocentric radial velocity of/? Pic is 20.0+0.5 km s -1 . 

7. The radial velocities of ions observed in absorption are 
close to or at the systematic velocity of /3Pic (to a few 
kms"'). 

8. The estimated radial density from Na I in emission predicts 
a column density similar to the one observed in absorption, 
meaning that most, perhaps all, Nal is distributed in the 
observed disk. We have no reason to doubt that originators 
of other 'stable' gas absorption components also belong to 
this extended disk. 

9. Our disk models show that assuming hydrogen to be the 
braking agent for metals pushed out by radiation pressure 
in the f3 Pic disk leads to contradictions with observations. 

A more detailed study of the chemical composition of the [3 Pic 
disk will be presented in a forthcoming paper. 
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